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Abstract 

The distinction between type I and type II superconductivity is re- 
examined in the context of the SO (5) model recently put forth by Zhang. 
Whereas in conventional superconductivity only one parameter (the Ginz- 
burg-Landau parameter k) characterizes the model, in the SO(5) model 
there are two essential parameters. These can be chosen to be n and 
another parameter, j3, related to the doping. There is a more complicated 
relation between k and the behavior of a superconductor in a magnetic 
field. In particular, one can find type I superconductivity, even when k is 
large, for appropriate values of (3. 



1 Introduction 

Several years ago, the SO (5) model was proposed as a description of high- 
temperature superconductors Q. At low temperatures these materials exhibit 
superconductivity (SC) or antiferromagnetism (AF) depending on the doping; 
the SO (5) model attempts to unify these two phenomena (both of which involve 
spontaneous symmetry breaking) into a single symmetry group. 

Not surprisingly, the SO (5) model is somewhat more complicated than the 
corresponding model describing conventional superconductivity (which is based 
on the smaller group SO (2)); correspondingly, richer behavior can be seen. In 
conventional superconductivity, one parameter, the Ginzburg-Landau (GL) pa- 
rameter a, determines the behavior. For instance, n = 1/V2 = n c marks the 
boundary between type I and type II superconductors. 

In the SO(5) model, there are two dimcnsionless parameters, which can be 
taken to be k and a second parameter, (3 (to be defined below); the latter is re- 
lated to the degree to which SO (5) is explicitly broken (by doping, for instance). 
As will be shown below, this second parameter has a profound influence on the 
magnetic behavior of an SO (5) superconductor; for instance, the critical value 
of k in the SO (5) model is given by k c ((3) = (1/^)^(1 + /3)/(l - (3). One sees 
that as (3 — > 1 _ , k c — > oo, which is a dramatic departure from the conventional 
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value. This is of some significance since high-temperature superconductors typ- 
ically have k f, and therefore are normally thought to be extreme type II 
superconductors. Such a conclusion would be premature in the SO(5) model, 
however, since (as the above relation shows) one could have k c >k»1. 

In this paper we will first remind the reader of the magnetic properties of 
conventional superconductors. The standard approach based on the surface 
energy density of a boundary between a normal phase at critical magnetic field 
H c and a superconducting phase will be briefly reviewed. We will then discuss an 
alternative approach based on vortex energetics. Although much of the above is 
fairly familiar material, this discussion will establish notation and set the stage 
for the parallel discussion in the context of the SO (5) model. 

Next, we will discuss the case of SO (5) superconductivity. After a brief 
review of the model itself, the critical fields H c and H C 2 will be calculated. 
From these, the critical value of the parameter n can be calculated explicitly as a 
function of (3. This expression can be confirmed numerically, either by examining 
the surface energy density or vortex energetics. In the latter approach, one can 
see that the dramatically different behavior in the SO (5) model is due to the 
possibility of an antiferromagnetic core in the vortex (2[ [| [| . 



2 Review of conventional superconductivity^] 

A conventional superconductor is described by the following Helmholtz free 
energy in the GL theory: 
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where /„ is a constant, h = V x A is the microscopic magnetic field and ai, b 
are parameters. The minimum of the potential is \4>\ 2 = af/b = v 2 . 

The Helmholtz free energy is minimized when the temperature T and interior 
magnetic induction B are constant (i.e., thermodynamic equilibrium). One uses 
the Gibbs free energy for equilibrium where T and the external magnetic field 
H are constant. These two quantities are related by the following Legendre 
transformation: 

g(H,T) = f(h,T)-±-hH. (2) 

The second term in this equation is the microscopic version of {Air) ~ l BH, since 
the magnetic induction B(r) = (h(r)). 

The critical magnetic field H c is defined as the field which establishes the 
following condition on the Gibbs energies in the normal and superconductor 
regions: 

g n (H c ,T) = g s (H c ,T). (3) 



1 For a more complete discussion on conventional SC, see |5| 
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We also have the following relations, the first one due to the Meissner effect: 

g s (H,T) = f s (T), 
H 2 

9n{H,T) =/„-—. 
From Eq. (If]) , we can obtain a correspondence with the parameters of the theory: 



/ n (T)-/ s (T) 
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One finds: 
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where $o = hc/e* is the quantum of magnetic flux, A = (m*c 2 /Ane* v 2 ) 1 / 2 is 
the penetration depth and £ = (h 2 /m*af ) 1//2 is the coherence length. 

There is also another critical field, H C 2, which is the strongest field allowed 
in a superconducting region before the system goes to the normal state. It is 
given by: 



H C 2 = 



$0 = V2kH c , 



2tt£ 

where, as usual, k — A/£. We see that if n is equal to its critical value l/\/2, 
H c = H C 2 (more on this in Sections 2.1 and 3.2). 



2.1 Surface Energy 

As is well known, a superconductor will behave in one of two ways when placed 
in a magnetic field. In a type I SC, macroscopic normal regions where h = H c 
will form, while in a type II SC, a lattice of vortices of flux $o will appear. 
One way to determine under which circumstances a SC is of type I or II is 
by studying the energy per unit area of a boundary between a normal region 
and a SC region. If this energy is positive, the system will tend to minimize 
the interface's surface area in order to lower this energy as much as possible, 
indicative of a type I SC. Similarly for negative surface energy, this surface area 
will be maximized, inducing type II behavior. 
Analytically, the surface energy is given by: 

$W T _H 2 C \ 
s 32tt 3 A 3 4tt ' 

where: 

1 = J°° ds\^-{\- f) + h 2 -V2hy (5) 
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Consider therefore a field configuration between the normal state for x — ► 
— oo and the SC state for x — > +00: 

tf> = (f>(x), A = A(x)y. 

Let us go to a description in terms of dimensionless variables, setting: 

ch - ch 
x = As, = vf, A = —a, h = —rh. (6) 

We obtain the following dimensionless free energy: 

F = J da {/„ + ±f' 2 + (a 2 + f l)f + h^. (7) 

Minimizing F yields the following equations (the so-called ID GL equations): 



7^ + (!-a 2 )/-/ 3 = 0, (8) 
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k 2 ds 2 

d 2 a 



ds2 -«/ 2 = 0. (9) 

The sign of the surface energy can be determined qualitatively in the limit 
of either small or large value of k in Eq. (||) , by looking at the length scales of h 
and /. By construction (see Eq. (||)), the length scale of h is lh = 1, while that 
of / turns out to be If = 1/(^/2k)). 



(a) Small k 



h 




Figure 1: Field solution behaviors for the conventional superconductor (a) type 
I (b) type II 

From Fig. |l|, we can see how the sign of the surface energy is determined 
by the behavior of / and h in those limiting cases. In the first limiting case 
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(Fig. ^]a) which corresponds to a type I sample (k -C 1, so that If ^> 1^), the 
energy comes from the expulsion of the magnetic field within the SC region. 
More explicitly, h is rapidly expulsed at the interface while / slowly assumes 
its asymptotic value. This gives rise to a positive scalar field energy density 
without a cancelling negative magnetic energy; the surface energy is positive. 
Fig. 0b shows the reverse case of a type II sample, where / quickly reaches its 
asymptotic behavior while the magnetic field slowly decreases, resulting in a 
negative surface energy. 

The critical value of k is the one for which the surface energy (or equivalently, 
the integral (||)) is zero. The well-known result n c = l/\/2 can be determined 
either numerically 0, or by making the oracular observation that when / and h 
obey the equations of motion for K = 1/ \/2, then the integrand of @ is zercfj 
This observation can be explained by a recasting of the work of Bogomol'nyi 
[|| (to be discussed in the next section) into a form appropriate to the surface 
energy. 



2.2 Vortex Energetics 

An alternative approach to determining the behavior of a superconductor when 
placed in a magnetic field is the study of the energy of vortices as a function 
of their winding number to. In this section, we will see what information can 
be drawn from the analytic expression, as well as an important limiting case, 
to — > oo. In this limit, we will also show how one can establish a link between 
the vortex and surface energies. 

As mentioned above, when placed in a magnetic field, the system will even- 
tually find the configuration which minimizes the free energy. We can construct 
a function J-(m) = F(m)/m which is the free energy of a vortex of winding 
number to divided by the winding number, or in other words the free energy 
per flux quantum of an TO-vortex. For large to, if !F(rri) is an increasing func- 
tion of to, it is more energetically expensive for a magnetic field to penetrate in 
vortices of large winding number; rather, a lattice of vortices of winding number 
1 will form. This, clearly, is what we expect of a type II superconductor. Sim- 
ilarly, if J-(m) is a decreasing function of to, the least energetic configuration 
will be a macroscopic normal region containing a large number of flux quanta; 
this is what we expect of a type I superconductor. 

For the to- vortex configuration, we use an ansatz that is generalized to wind- 
ing number to: 



e* r 

Again we can use the dimensionless variables defined in (^). Here, however, s 
denotes the radial variable in a 2D plane, as opposed to the cartesian coordinate 

2 A numerical calculation of the surface energy is eiven in H. 

3 Note that there is no need for the integrand of (H) itself to be zero; the weaker condition 
that the integral is zero would suffice. Indeed, we willsee in the next section that in the SO (5) 
model the analogous integrand is not zero when the integrated surface energy is zero. 
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orthogonal to the domain wall. 

Using these changes in Eq. ([[]) yields the vortex energy: 

,.. J f*f{i + ( rf + £) , + ^( / - + (= + »)V 

-(i-/ 1 )/ 2 



as well as the equations of motion: 



K z \ S \ S 



/ + / - / 3 = 0, 



h'+(—+a) f = 0. 

\KS I 

It is possible to rewrite Eq. ( |io|) in the following form 



(10) 
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where x = ko(/ 2 — 1) +mf 2 /s. Integrating the last two terms yields the winding 
number to, independent of the details of the field. Clearly, if k = 1/V2, the 
third term vanishes and F consists of the winding number to plus two positive 
semidefinite terms; it will be minimized if they are zero. One can see that these 
terms are in fact zero if / and a satisfy the equations of motion. In this case, 
T{m) = F(m)/m = 1 independent of m, which indeed corresponds to critical 
behavior with regard to vortice lattice stability, as there is no preference for 
either type I (coalescing) or type II (lattice) behavior. 

We can determine the behavior of T(m) numerically for a large value of the 
winding number to. Our numerical results already suggest that J-(m) reaches 
asymptotically a constant value at to — > oo. In this limiting case, the core of the 
vortex is very wide, and the fields fluctuate rapidly at only one spatial position- 
that is at the vortex perimeter. Let us define this position as s m . Analytically, 
we will see this more clearly by inserting the limiting forms of the field into 
Eq. (HI), which are: 



s <C s ni 

S > Sr 



/->o, 

/->!, 



~s, h = 0, 

to , a 1 



After integration, the second and third terms of Eq. (J13|) become, respectively: 



kH, 
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We have to and the remaining bracketed term, which does not contribute 
to the integral (this, of course, stems from our choice of step-functions as the 
role of fields / and h, since we essentially remove the presence of the surface 
altogether). Let us rewrite the variable s m in terms of the quantum of magnetic 
flux $0; using our m — > oo ansatze: 



Therefore: 



to J to 



, 2m 

si - 



kH c 



We find: 



Tim) = _ (to — > oo). 

If we add thickness to the interface, the first term of Eq. (E3h, which was 
vanishing, will instead be proportional to yfm (i.e., a "perimeter" term cx s m ). 

We are now able to establish a correspondence between the free energies of 
the vortex and wall domain configurations. Indeed, for a large value of to, the 
perimeter of the vortex is large enough to significantly lower the curvature of the 
vortex locally, likening a given small region into the form of a domain wall. We 
should expect to see confirmation of this, using the appropriate approximations. 
We will do this here. 

Let us start from the Gibbs energy of the vortex and derive an expression 
that clearly shows the relation with the domain wall Gibbs energy. We have: 

Gv =F V - g$ oTO = ( 2k 2 F(to) - VW) . 

Substituting F(m) with Eq. (|l3| ) and expanding the terms in square parentheses: 

G v = 7^— |tj (m+ r sds (^(1 - J 2 ) 2 + Kh(l f) 



16tt 2 A 2 V Jo { 2 
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s / \ s 



f' 2 -2ff l- + na + - + na) f - V2 K m . (14) 



We notice that the following terms are equal to a familiar quantity: 



s I ds 

which, once integrated, cancels the to term. If we take Eq. ( |TT| ) and multiply it 
by /k 2 , we have: 

/ - + na) 2 f 2 = ff' + K 2 (f-f). (15) 
s / 
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Next, we change variables s — ► z = s — s m , where it is understood that the 
region of interest is \z\ <C \s m \. Therefore 

/OO pOO 
(z + s m ) dz rj / s m dz. 
- oo J — oc 

Substituting ( p"5| ) into ( |l4| ) and simplifying, we have: 



167 



.2 



dz 



dz . (16) 



We recognize the first integral: it is s m times I, as defined by Eq. (||); 
the second integral is zero due to the boundary conditions. Therefore, the link 
between the Gibbs energies between the vortex and domain wall case is summed 
up by the following relation, valid for large m: 



— tt, — 9,9-^ = G s (2ir\s m ). 



In this section, we have studied vortex energetics as an alternative approach 
to the determination of the type (I vs. II) of a given superconductor, in the 
conventional case. We now turn our attention to the case of SO (5) supercon- 
ductivity. 



3 SO (5) Superconductivity 

Having briefly reviewed the basics of conventional SC, let us turn to the case of 
SO (5) superconductivity. As mentioned in the introduction, the SO (5) model 
attempts to unify the superconducting and antifcrromagnetic aspects of high- 
temperature superconductors. The order parameters for these two phenom- 
ena form a 5-dimensional real field whose dynamics has an approximate SO (5) 
symmetry, (fa, fa, r)i, r] 2 , where <p = fa + ifa and r\ = (rji, r\i, r] 3 ) are, re- 
spectively, the order parameters of superconductivity and antiferromagnetism. 
This so-called superspin is described by a free energy which will be given be- 
low; explicit SO (5) breaking enables spontaneous breaking of either the SO (2) 
symmetry of <j> or the SO(3) symmetry of 77, describing SC and AF, respectively. 

Mathematically, the low-energy effective theory is described in terms of the 
following free energy jjj : 



F = I dx I ^~ 



8tt 2m* 
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(17) 



where h = V x A is again the microscopic magnetic field. 
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The nature of the ground state can be determined by examining the poten- 
tial. We consider the most general symmetry-breaking potential including even 
powers of the helds up to quartic terms, 

a\ 2 a 2 2 b 1 4 + 2b 3 2 ? ? 2 +b 2 7 ? 4 
V(<t>,rj)=-—cl) - —t] + , 

where we have written <fi = \<f>\ and rj = |ry|. This potential is invariant under an 
5*0(3) x 5*0(2) symmetry. S0(5) symmetry would be attained by setting the 
two mass parameters and the three quartic coupling constants to be equal; for 
an approximate symmetry their values are approximately equal. To simplify the 
analysis we will set the quartic parameters to a common value b\ = 62 = &3 = b; 
explicit symmetry breaking in the potential is found only in the quadratic terms. 

Since we are interested in the magnetic properties of the superconducting 
state, we restrict ourselves to parameters which describe a SC ground state. 
This is done by requiring that the global minimum of V be on the <f> axis. The 
ground state will then have the value (<f>, 77) = (v, 0), where v = ai/y/b; the 
global condition is fulfilled if (3 = a^/a 2 < 1. Note that we recover SO (5) 
symmetry if (3 = 1 and, therefore, this is the value corresponding to critical 
doping. For (3 > 1 the ground state is antiferromagnetic. 

The goal of this section is to calculate the surface energy, whose sign indicates 
the type of superconductor described by the model. We will first set up the 
physical context of the surface energy and derive the critical fields H c and H C 2- 
These enable us to determine the value of function of (3. Subsequently, 

we will give numerical calculations which confirm our analytical results. 

3.1 One Dimensional Free Energy 



As outlined in Section 2.1, the surface energy is the Gibbs free energy per unit 
area of a domain wall separating a normal (non-superconducting) phase in the 
thermodynamic critical field H c and a superconducting phase in the absence 
of a magnetic field. We use the same ansatz for </> and A as mentioned in the 



beginning of Section 2.1, and with the AF order parameter written 77 = rj(x)x, 
the free energy ( |l7| ) becomes: 



F = / dx lf n + — 




(18) 



It will again be useful to use the dimensionless quantities defined in (|E|) . Writing 
r)(x) = vn(s), we find, after some algebra, the following dimensionless free 
energy: 

F = f dx {/„ + ^(/' 2 + n' 2 ) + a 2 f + h 2 - f - (3n 2 
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~(/ 4 + 2/V+n 4 )}. (19) 



Minimizing with respect to /, n and a, we obtain the SO (5) GL equations: 

f-f- n 2 f = 0, (20) 



111 

k 2 ds 2 



a 2 ) 



ds 2 



af, (21) 



1 d 2 n , 9 , 

^d^ +(3n ~ n Uf °- (22) 
Note that if we set n — in the first two equations, we recover the conventional 
GL Eqs. as expected. Note also that we are now left with two parameters 

(k, (3) rather than three (a 2 , a 2 , and b). The first of these two is the traditional 
GL parameter k = A/£, which is usually quite large for high T c SC (typically 
between 15 and 100). The second fundamental parameter is (3, which is related 
to the doping of the system away from critical doping (that corresponding to 
the AF/SC transition). /3 can be written as a function of the chemical potential 
/z as P = 1 — 8m*£ 2 x(/i4 2 — ^t 2 )/?i 2 , where x is the charge susceptibility and /j c 
is the critical value of the chemical potential Q . 

3.2 Calculation of H c and H c2 

When a critical magnetic field is applied in a region, the SC order parameter will 
be forced to zero, destroying superconductivity. In order to minimize the Gibbs 
free energy, a nonzero AF order parameter will be induced. This is because, / 
being held to zero by the magnetic field, the potential in this region is minimized 
if n = y/]3, as can be seen from the potential terms in (|l^). The value of the 
critical field is determined by the competition between the positive potential 
energy of the AF state (relative to the SC ground state), on the one hand, and 
the negative magnetic energy (—H 2 /8n), on the other. 

Let us compute, first of all, the critical magnetic field H c . The result differs 
from the conventional one because an AF order parameter will now appear in 
the normal region. 

Far into the superconducting region (x — > oo), we get the following behavior 
for the fields: h = 0, / = 1, n = 0. Thus from Eqs. (@) and @, 

g s (H,T) « g(h = 0,T) = f s (T) = f n (T) - ^. (23) 

In the normal region (x — > — oo), we include a nonzero AF order parameter, and 
set h = H, n = v 7 /?, / = 0, giving, 

9n (H,T) = f n (T)-f b -^. (24) 
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By equating the two Gibbs free energies, we find the following /3-dependent 
expression for the critical magnetic field: 



H C (P) = H° c ^/T-~fP, (25) 

where H® is the conventional expression (Q). Note that H c {(3) will refer to the 
critical field in the SO (5) model, and must be distinguished from the one in 
the GL theory (from now on, we will write H®). Eq. ( ^5|) tells us that for 
(3=1, H c {(3) vanishes. In fact, this value of (3 corresponds to the critical value, 
where the transition from SC ((3 < 1) to AF ((3 > 1) occurs; it is natural that 
any magnetic field, however small, will break SC, since SC is not energetically 
favored over AF: their free energies are equal. As (3 decreases, H c (f3) reaches 
a maximum when (3 = 0, i.e., when the sample is at maximal doping. In this 
case, the value of H c ((3) is the same as in conventional GL theory with the 
same n. Thus, far in the normal state, the critical field will take the value 
H C (J3) = (l/v^V 1 - /3 2 since H° = 1/V2 in the conventional SO(2) case (see 
Section 

Also of interest is the critical field H C 2- For a type I SC, at T < T C) when 
H > H®, the SC phase is destroyed. After that, if H is reduced, the SC is 
recovered at H C 2 < H° (phenomenon of supercooling). For a type II SC, H C 2 
corresponds to the maximum magnetic field in which a vortex lattice can appear; 
it is always higher than H®. Then, H C 2 is the strongest field which permits the 
SC phase to occur. 

Near H c2 , the order parameter is small (the system is near the transition), 
and we can use the linearized SO(5) GL equations. On the other hand, we 
can show that, at first order, the microscopic field h can be taken equal to the 
applied field H. Beginning with the free energy Eq. (|l7|), we choose the Landau 
gauge H = Hz and Ay = Hx. Linearizing around <j> = and 77 = a^/^/b (the 
values of fields at the transition), the equation for the SC order parameter takes 
the form Hi: 



-2 



—Hx— ( 27T ^ X \ 
$0 X dy ^ $ j 



^ = -^(l-/3)V, 



where we have set H = X 2 H , <fi = vip and $0 is the flux quantum. 

This equation can be put in an harmonic oscillator equation form for x, by 
scting ip = f(x) exp (—i(k y y + k z z)) and the quantization of the energy leads to 
a set of allowed values of the magnetic field. The highest possible value of H n , 
then, corresponds to the upper critical field and we find 

H C 2(f3) = ^(l-l3) = V2 K H° c (l-f3). (26) 

Note that it is the conventional critical field that appears in the right hand of 
Eq. ©. 

In the GL theory of conventional SC, we have seen that the boundary be- 
tween type I and type II SC occurs when k = n c — l/v2; at this value, 
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H C 2 = H®. We can already anticipate the value of k c in the S0(5) model 
as the value for which these two critical fields are equal. 
From ( p5| ) and (pfj) we find: 

k c ((3) = (l/V2)v/(l+/5)/(l-/3). (27) 

Thus we are able to obtain an analytical expression for function of the 

doping. Two facts emerge from Eq. (p7|): firstly, we recover (as we must) the 
conventional result when (3 = 0, and secondly, for any f3, a value of n c exists 
which separates type I and type II behavior (Figure |2|). 

As (3 approaches 1 (the SO (5) symmetric case), the instability of a vortex 
lattice (in the type I region) will be amplified, and therefore, to reach the type 
II region, we will need a higher k. In practice, it is easier to vary the doping [3 
(or n) than the GL parameter k in a sample. 

We can reverse this point of view: if we consider a SC at a fixed value of n, 
we will find that, for (3 > (3 C (defined implicitly in (]27])), the sample will exhibit 
type I behavior. Thus, for high values of n (i.e., n ^> 1), we also need high [3 
(i.e., [3 approaching 1) in order to see type I behavior. We will see in the next 
section that the numerical calculation of the surface energy agrees perfectly with 
these results. 



Limit between type I and type II 
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Figure 2: k c as a function of the doping f3. The curve delimits the different ty pes 
of SC. The points come from numerical calculations (see the end of Section [^3] 
and Table 1). 



3.3 Surface Energy 

In this section, we will derive an expression for the surface energy in the SO (5) 
model, and rederive from it numerically the function n c {(3)- We will begin by 
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analyzing qualitatively the surface energy; we will sec that for k, sufficiently large 
(small), the surface energy is negative (positive). However, even this qualitative 
argument indicates that the critical value of k depends on (3, as anticipated in 
the previous section. We will then present numerical results since, as in the 
conventional model, the surface energy cannot be evaluated analytically; these 
numerical results agree with the conclusion of the previous section based on the 
critical fields. 

To analyze qualitatively the surface energy, consider Fig. |^. In the normal 
region as x — > — oo, the positive contribution to the energy coming from the 
scalar fields / and n exactly cancels the negative condensation energy of the 
magnetic field. In the SC region, as x — ► oo, the scalar field and magnetic 
energy densities individually are zero. Inside the surface, variations in the fields 
contribute to the surface energy. The details of these variations will determine 
the sign of the surface energy. 

We can characterize the fields' variations by their characteristic lengths. 
From Eqs. (pofj^ ), these are seen to be 

lh = h l f = -/7T> ln = — r W ( 28 ) 

V2k - [3 

At issue is whether in going from one asymptotic regime to the other, the 
scalar fields or the magnetic field are capable of changing more quickly. If the 
scalar fields vary faster, the magnetic field will penetrate into the superconduct- 
ing region, and the negative magnetic energy density in this transition zone will 
not be accompanied by a positive scalar energy density, resulting in a negative 
surface energy. If the magnetic field varies faster, the situation is reversed and 
the surface energy will be positive. Thus we need to compare the larger of l n 
and If with fo. We see from Eq. (^8j) that l n is always greater than If, so l n 
determines how quickly the scalar fields can vary. Two cases are important to us 
when k l/v2- First, when (3 is small (Fig. ||a), we find from ( p8[ ) that lh > l n , 
and as described above, the surface energy will be negative. This corresponds 
to the behavior of a type II SC. 

On the other hand, when (3 is large (Fig. ||b), we have the opposite situation: 
lh < l n , and the surface energy will be positive. This is the case of a type I SC. 

By continuity, for any k, there must exist a (3 where the surface energy is 
zero. This provides an alternative way of finding ft(/3), and must agree with 
Eq. ©. 



We define the SO(5) surface energy in the same way as in Section 2.1. Sub 



tracting Eq. (24) from Eq. (H) we are led to 



* / H'--^-U»-3-^]>- f2!,i 




Using (hq) we find, after some algebra, 



0\2 
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(a) p small 



(b) p large 



Figure 3: Field profiles in surface for k » l/v2. (a) (3 small (b) (3 large, 
where 



1 = j°° ds - / 4 - n 4 ) - /V + h 2 - V2h^~W< 



(30) 



with A and H® defined in Section |^. Therefore, the (3 dependence appears only 
in /, and the expression reduces to the conventional one, Eq. (|5|), when we set 
n = = 0. 

Of interest are the conditions under which this integral is zero. It will be 
recalled that in the SO(2) case, when k = l/\/2, not only is the integral zero 
but the integrand is zero - a stronger statement, and one which is explained by 
Bogomol'nyi One might hope that the same feature would be found in the 
SO (5) model, along with an analogous analytical explanation. Unfortunately, 
such is not the case, and the surface energy must be studied numerically to find 

We have solved Eqs. (|2C|-p2[) numerically, using a relaxation algorithm. The 
boundary conditions appropriate for the problem are: 



x^oo: /->1, h^O, n^O. lj 

Having found the solution, it is easy to compute the dimcnsionless surface energy 
/ as a function of the doping j3 by numerical integration of (^). We present, in 
Fig. ||, the results for k — 0.707 and k = 7.07 respectively. The case k < I/a/2 
is of no interest since the surface energy is positive for all (3. 

Fig. |^a corresponds to n c for conventional SC. First of all, we see that the 
surface energy is zero for f3 = 0, as expected since, at this particular value, we 
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recover the SO(2) case. For all other values of 0, we have a positive surface 
energy. Thus, the sample remains in the type I region for any value of the 
doping. This is in agreement with Eq. (p7j), as seen in Fig. |^. Note also that 
the energy reaches a maximum and decreases until it reaches zero again at 
0=1. At this point, h = throughout the sample (see Eq. (|3l|)), and the 
ground state is SO (5) symmetric. Hence, no broken symmetry appears and for 
topological reasons (since the surface energy can be treated like a soliton), the 
surface energy must be zero. 
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Figure 4: Surface energy 1(0) for a) k = 0.707 and b) k = 7.07. 

A more interesting case is presented in Fig. |]b, corresponding to k = 7.07 
(in the conventional type II superconductor). The curve starts at — 

at the conventional value and passes through zero at sa 0.98, in agreement 
with the analytical result (|27|). Below this critical value, the sample has type 
II behavior since the surface energy is negative. Above, the surface energy is 
positive, corresponding to type I behavior, i.e., vortices in this region are not 
energetically stable. Observe that the surface energy is again zero for = 1. 

We can determine in this way n c as a function of 0; the results are reported 
in Table 1, and are in excellent agreement with Eq. (p7|). 



a) k = 0.707 
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0.3945 
0.5288 
0.6620 
0.7887 
0.8976 
0.9410 
0.98 
0.9998 



k c ((3) (numerical) 



1 

1.2 
1.5 

2 

3 

4 
7.07 
70.7 



From (|27|) 



1.07 
1.27 
1.57 
2.06 
3.04 
4.06 
7.04 
70.71 



Table 1 : Comparison of k c obtained numerically via vortex energetics and 
surface energy, and analytically via Eq. ( ^7|) . 



3.4 Vortex Energetics in the SO(5) Model 

As discussed in Section ^j^, the magnetic behavior of a superconductor can also 
be determined by studying vortex energetics [Q. Thus, one can calculate the 
free energy per quantum of magnetic flux of an m- vortex J-(m) = F(m)/m; if 
this function is of positive or negative slope, the superconductor is of type II 
or I, respectively; zero slope indicates the critical case, which can be expressed 
in terms of k c ((3). In the SO (5) model, the numerical calculation of jF(m) for 
n = 7.07 and 70.7 are given, respectively, in Fig. 2c and Fig. 2d in From 
the first of these, we can see that for (3 < 0.98 Tim) is of positive slope, while 
for (3 > 0.98 it is of negative slope; to a very good approximation, the slope 
is zero at (3 = 0.98. Thus, k c {(3 = 0.98) = 7.07. Similarly, Fig. 2d of the 
same reference indicates n c ((3 = 0.9998) = 70.7. A c omparison of this method 
with the analytical method discussed in Section 3.2 is shown in Table 1; the 
agreement is excellent. 



4 Summary and Conclusions 

In this paper, magnetic properties of superconductors were discussed, with an 
eye towards analyzing the SO (5) model of high-temperature superconductivity. 
The Ginzburg-Landau theory of conventional superconductivity was discussed in 
Section ^ mainly to review some fundamental and well-known results and also to 
establish notation. The energy density of a surface separating superconducting 
and normal regions at critical applied field, and its importance for determining 
the behavior of the superconductor, were discussed. In addition, we reviewed 
the approach of Bogomol'nyi wherein the vortex energy is written in a form 
that shows clearly that at k c , the energy is proportional to the winding number 
of the vortex. Finally, we discussed an alternative approach to determining the 
critical field which is suggested by Bogomol'nyi's work, namely, studying the 
energy of an m-vortex as a function of its winding number. Both analytically 
and numerically, this approach can be shown to be equivalent to the surface 
energy approach. 
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In the second part, we introduced briefly the SO(5) model, and derived, in 
the context of that model, expressions for the critical magnetic fields, H c ({3) and 
H C 2(P)- These enable us to calculate an analytic expression for k c ((3). Next, 
we studied the surface energy of a boundary between superconducting and non- 
superconducting regions in an applied critical field. Finally, we discussed briefly 
an alternative approach based on vortex energetics. 

It is worth comparing the two numerical approaches to calculating k c , namely 
by calculating the surface energy or by vortex energetics. In the former ap- 
proach, one must solve the field equations (for instance, (|2fj|-p2|) in the SO (5) 
model), evaluate the surface integral (|3^), and vary [3 for fixed k (or vice- versa) 
until the surface integral vanishes. In the latter approach, one must solve the 



vortex field equations (the generalization of fllip and ( |12| ) to SO (5)), evaluate 
the energy as a function of winding number, and determine for which k, f3 the 
function T{jn) is independent of m for large m. The latter method is somewhat 
more demanding numerically, essentially since several values of m must be con- 
sidered for each k, (3, though it is perhaps more intuitive, since one can think 
of type I behavior as due to an energetic preference for the coalescence of many 
single vortices into one large m-vortex. 



The main new result in this paper is Eq. (27), which expresses the value of 
k as a function of (3 corresponding to the boundary between type I and type II 
superconductors. This result is rather important: it states that as j3 — ► 1 (i.e., 
as the doping is reduced to its critical value), then even with a large value of 
k, the material would be a type I superconductor, contrary to what one would 
naively expect. Some implications of this have been discussed in [Q. 
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